The ground state of a mixture of two species of fermionic atoms in ID optical lattice 
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In this paper, we investigate the ground state properties of a mixture of two species of fermionic 
atoms in one-dimensional optical lattice, as described by the asymmetric Hubbard model. The 
quantum phase transition from density wave to phase separation is investigated by studying both 
the corresponding charge order parameter and quantum entanglement. A rigorous proof that even 
for the single hole doping case, the density wave is unstable to the phase separation in the infinite U 
limit, is given. Therefore, our results are quite instructive for both on-going experiments on strongly 
correlated cold-atomic systems and traditional heavy fermion systems. 
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I. INTRODUCTION 

Rapid progress in Bose-Einstein condensates in optical 
lattices [l|, [2|, y, U, [5| has opened fascinating experimental 
possibilities in condensed matter physics, atomic physics, 
and quantum information. For example, the experiment 
on neutral atoms trapped in the periodic potential of 
an optical lattice has been used to realize an array of 
quantum gates Moreover, cold-atomic systems are 
intrinsically related to many-body models in condensed 
matter physics. Compared to solid state systems, cold- 
atomic systems could be better experimentally control- 
lable. Thus, the investigation of cold-atomic systems will 
not only help us to have a deep understanding of known 
physical phenomena in many-body systems, but also pro- 
vide hints to explore new areas of physics. Such a beau- 
tiful prospect has attracted many theoretical attentions 

UMMMMM- 

Quite recently, experiments on fermionic atoms 
trapped in optical lattices Q were carried out which 
opened a door for us to find deeper insights into some 
essential problems in condensed matter physics, such as 
BEC-BCS crossover, superfluidity, and Mott-insulating 
phase. It was proposed that ultracold fermionic atoms 
exposed to the periodical potential of an optical lattice 
could be an ideal realization of the Bosc- Hubbard model 
[f|, the spin-dependent Hubbard model @, and the an- 
tiferromagnetic states or d-wave pairing states 0. The 
unique control over all relevant parameters in these sys- 
tems [5] allows people to carry out experiments which 
are not handy with solid state systems, so they marked a 
milestone towards the understanding of some fundamen- 
tal concepts in quantum many-body systems. 

In this paper, we consider a system of two species 
of fermionic atoms @ with equal numbers (or one type 
of fermionic atoms with spin-depent hoping integral [7() 
away from half-filling in an one-dimensional optical lat- 
tice, as described by the asymmetric Hubbard model. 
0, EDI G3 (AHM). The system is expected to have a 
density wave (DW) state and phase separation (PS) of 
two atom species state [l5[ , and we investigate the quan- 
tum phase transition (QPT) from the DW state to the 



PS state in this system by studying both the quantum 
entanglement and traditional DW order parameter. We 
show that the entanglement can help us to witness crit- 
ical phenomenon and shows scaling behavior around the 
critical point. The phase transition is also clarified by 
the competition between two different modes of struc- 
ture factor. A global phase diagram as a function of the 
local interaction U and the ratio of two hoping integrals 
is then obtained under different filling conditions. More- 
over, we give a rigorous proof, that even for the case of a 
single hole doping away from half-filling, the DW state is 
unstable to the PS state in the infinite U limit. As will 
be shown below, if we regard two regions in the PS phase 
as one solid-like region of heavy atoms and another as a 
liquid-like region of light atoms, respectively, the QPT 
is just a physical realization of the quantum solvation 
process [16| in the optical lattice. Therefore, our results 
are quite instructive for on-going experiments on strongly 
correlated cold-atomic systems. The behavior of entan- 
glement in this system can help people to have a deep 
understanding of the critical phenomenon. 

This paper is organized as follows. In section HH we in- 
troduce the Hamiltonian of the AHM, and show how to 
realize the model in the quasi one-dimensional periodical 
potential of an optical lattice. We will also briefly in- 
troduce the background of the model in the condensed 
matter physics. In section IIII1 we study the ground- 
state entanglement of the system. We will show that a 
schematic phase diagram can be obtained from the entan- 
glement between a local part and the rest of the system of 
a finite sample. In section ITVl by studying the structure 
factor of the density distribution of heavy atoms, we can 
obtain a quantitative phase diagram for different filling 
conditions via both the exact diagonalization (ED) and 
density matrix renormalization group (DMRG) methods. 
In section |V1 we will give a rigorous proof that even for 
the case of a single hole doping away from half-filling, 
the DW state is unstable to the PS state in the infinite 
U limit. While if U is very large, the critical point then 
is approached linearly with 1/U. In section IVTl we will 
discuss the mechanism of the existence of the PS, the 
possibility of the PS in high dimension, and conditions 
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for experimental realization. Finally, we summarize our 
results in section fyTTl 



II. THE MODEL HAMILTONIAN 



The one-dimensional AHM is defined as 



i=l <5=±1 



In Eq.fl]), t a {a — a, (3) distinguishes the species of 
fermionic atoms (e.g., 6 Li and 40 K), ct g and Cj j(n a =|, | 
are creation and annihilation operators for a atoms at 
site j respectively, and n a — cJ.c CT , while U denotes the 
strength of on-site interaction. In this model, the Hamil- 
tonian has U(1)®U(1) symmetry for the general t a , and 
the atoms number N a = . rij ta , Np — J2j n j,P are con ~ 
served respectively. The total number of atoms is given 
by N = N a + Np, and the filling factor is n = N/L. 

The asymmetric Hubbard model |T]) can be used as an 
effective model to describe a mixture of two species of 
fermionic atoms in an one-dimensional optical lattice. In 
order to have a quasi one-dimensional system, we suggest 
that the optical lattice potential takes the form of 

V(x,y,z) = Vbsin 2 (fcr) + V± [sin 2 (ky) + sin 2 (fcz)], 
H 2 k 2 



Vn = V- 



V ± =uj 



2m 

h 2 k 2 



2m 



(2) 



Here k = 2ir/A and A is the wavelength of the laser light, 
and Vo and V± denote the maximum potential depth 
along the x direction and in the yz plane respectively. 
The potential depth is measured in units of the recoil 
energy h 2 k 2 /2m. In order to freeze the hoping process 
in the yz plane we should have V± 3> Vq- For a single 
atom in the periodic lattice, its wave function is the Bloch 
state, which is actually a superposition of well localized 
Wannier state. Therefore, if we restrict ourself to a very 
low temperature, where the thermal fluctuation cannot 
excite the atom to the second band, the Wannier state 
can be approximated by the ground state of a single atom 
in the potential well. For the present case, the ground 
state can be written as 

* (x, y ,z)^(^) 1/2 e-^^^ (,), (3) 

where 

fmu)\ 1 / i j 
*o(*) = (— ) e *• 

hk 2 y/UT hk 2 ^ 

u± = . u; = . (4) 

m m 

Then the hoping matrix element between the two adja- 
cent sites i, j can be calculated as 



t = - J drw a (r - ri) 



ft 2 V 2 
2m 



+ V]V>„(T-T j ), (5) 



which results in the hoping integral along the x direction 

h 2 k 2 



t : 



2m 



(6) 



Moreover, if two atoms, a and (3 occupy the same site, 
they will repel each other. The on-site interaction can be 
approximated with 



U 



4?rh 2 a 
^m a m (i 



K(r)| 2 K(r)| 2 dr, 



where a is the scattering length. Using the wave function 
of Eq. we obtain 



U 



4nh 2 a kv^Pv 1 / 2 

^Jm a m(3 y/TT 7T 



(7) 



Finally, if we have a system of two species of polarized 
fermionic atoms in the optical lattice, the hoping integral 
and the on-site interaction will have the form (in units of 

ta) 



tp_ _ ma 

U \d>ayJ-Km a /mf3 v 1 / i v 1 / 2 

ta 



A 



(v 7 ^ + 2y^T) 



Taking Li(a) and K (/3), two species of atoms, as an 
example and v±_ = 16, we have 



tp ~ 0.15, 
jj „ 35.87 a 



w i/4 



A + 8) 



(9) 



In condensed matter physics, the asymmetric Hubbard 
model is one of the most simplest two band models which 
is believed to describe many essential physical properties 
of strongly correlated systems. To understand the in- 
teresting phenomena which may happen in the ground 
state of the Hamiltonian (fTJ), it is very useful to look 
into the two limiting cases of Eq. (TB). If t a = tp, the 
AHM becomes the Hubbard model In ID, the Hub- 
bard model can be solved exactly by the Bethe-ansatz 
method. The wave function and the energy spectra 
then can be calculated exactly. In the large U limit, 
the Hubbard model can be approximated by the famous 
t — J model, in which the spin-spin interaction is of the 
antifcrromagnetic type. Therefore, it is widely accepted 
that the ground state of the Hubbard model at half- 
filling shows the spin-density wave. On the other hand, 
if tp = 0, the AHM becomes Falicov-Kimball model. 
[3 Eq|, [Hi Hi| In ID, it has been pointed out that the 
system will segregate into an empty lattice (with no (3 
atoms and all a atoms) and a full lattice (with all 
atoms and no a atoms) in the large U limit when away 
from the half-filling. Therefore, the two limiting cases of 
the AHM belong to different universality classes, a phase 
transition from PS to DW is expected to appear some- 
where on the U — tp plane. 
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FIG. 1: (color online) The changes of symmetry in the ground state wave function is analyzed by considering the quantum 
correlation, i.e. entanglement, between local block and other parts of the system. Here L = 6, N a = Np = 2, 1 = 2, and the 
anti-periodic boundary conditions are assumed in order to avoid level-crossing in the ground state. Four figures correspond to 
different block size: &{l = 1), b(Z = 2), c(7 = 3), and d(l — 4). 



III. GROUND STATE ENTANGLEMENT 

In recent years, studies on the role of entanglement 
in the quantum critical behavior [23[ have established 
a bridge between quantum information theory and sta- 
tistical physics. [24|, |25|, [26| It is believed that the en- 
tanglement, as a kind of quantum correlation, can help 
us identify quantum phase transition in many-body sys- 
tems. To have an intuitive picture of the global phase di- 
agram, we must first compute the entanglement between 
a local block and the rest of the system. For the present 
model, the local states on each site have four possible 
configurations, denoted by 

<^ = |0)», |/3>, h9>; I = 1,2,3,4. 

The Hilbert space associated with the L-site system is 
spanned by 4 L basis vectors. If we choose the periodic 
boundary conditions for N = An + 2 and antiperiodic 
boundary conditions for N — 4n, where n is an integer, 
the ground state is nondegenerate. Considering the re- 
duced density matrix of a block of I successive sites of 



the ground state 

pi =tr r |*)<*|, (10) 

the von Neumann entropy E V {1), i.e. 

E v (l) = -tr[p i log 2 (p / )] (11) 

measures the entanglement between the I sites and the 
L — I sites of the system. Like the well known fact in clas- 
sical optics that the three-dimensional image of one ob- 
ject can be recovered from a small piece of holograph due 
to the interference pattern of the reflected light beams 
from it, quantum superposition principle also allows us 
to see a global picture of the system from its local part 
[25l |2o| . As was shown for some typical models in con- 
densed matter physics, such as the extended Hubbard 
model[26[, the entanglement of the ground state can give 
us a global view of the phase diagram. 

^From this point of view, we show a three-dimensional 
diagram and its contour map of the entanglement with 
block size I = 1, 2, 3, 4 for 6-site system with N a = Np = 
2 in Fig. [T] It has been pointed out that in the extended 
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FIG. 2: The entanglement (left) and its first derivative (right) 
as a function of tp for various block size and a specified U = 
200. The inset shows the scaling behavior of the minimum 
point of dE v (l) /dtp at the critical point. Here L = 10, JV a = 
Np = 4. 



Hubbard model, the single-site entanglement can distin- 
guish the three main phases in the ground state. The 
reason is that the density distributions of the different 
modes in the reduced density matrix of a single site, such 
as the double occupancy, in the extended Hubbard model 
are sensitive in the quantum phase transitions. However, 
from Fig. QJa), the single-site entanglement in AHM is 
rather trivial in the large U region. It is not difficult to 
understand this phenomena. For the present model, the 
reduced density matrix of a single site has a simple form, 

mm 

pi = z\0)(0\ + u+\a)(a\+u-\(3)((3\+w\a(3)(a(3\, (12) 



in which 



", and w are the density distributions 



for different local states, and can be calculated as 

w = (n a np) = tr(n a nppi), 
u + = (n a ) — to, u~ = (rip) — w, 
z = 1 — u + — u~ — w. 



(13) 



In the large U limit, the double occupancy of two atoms 
on a single site is forbidden, i.e. w ~ 0. Then for a finite 
system with periodic or anti-periodic boundary condi- 
tions, (n a ) and {rip) are constants. This fact leads to 
a constant single-site entanglement during the evolution 
of tp in the large U region (Fig. QIa)). For the case of 
n a = 1/3, rip = 1/3, E v (l) has the value log 2 3, so the 
single-site entanglement is insensitive to the phase tran- 
sition from DW to PS. This property is very similar to 
the single-site entanglement in spin models [2 71 ] and the 
ionic Hubbard model. [28j 

Clearly, the transition from DW to PS is intrinsically 
related to the change of density distribution of one species 
of atoms on the lattice. In order to contain enough infor- 
mation of the density-density correlation from the point 
of view of the entanglement, more sites should be in- 
cluded into the block. According to this point, we show 



the two-site entanglement as a function of tp and U in 
Fig. [1] (b) , from which we immediately notice two differ- 
ent regions: one is an altiplano marked with warm color 
(denoted by "PS" in the contour map of Fig. [U(b)), while 
the other is a plain with cold color (denoted by "DW" 
in contour map of Fig. [T] (b). Taking into account the 
known fact of the two limiting cases of this model, such 
an obvious difference witnesses the critical phenomenon 
between two universal classes. 

In order to understand this obvious difference of the 
two-site entanglement in two phases, let us have a look 
at the structure of the corresponding reduced density ma- 
trix. For the AHM, the total numbers of a atoms and 
(3 atoms are good quantum numbers, which leads to the 
fact that for arbitrary block size I, there is no coherent 
superposition of local states with different values of N a 
and Np in the reduced density matrix. That is, the re- 
duced density matrix must have the block-diagonal form 
classified by both N a (l) and Np(l), i.e. 

Pl = diag{ P/ (0, 0), Pl (l,0), p,(0, 1), . . . , pi(l, I)} (14) 

where pi{n a ,rib) is a matrix which has n a a atoms and 
rib (3 atoms. According to the definition of von Neumann 
entropy (hence the entanglement), its magnitude is really 
determined by the distribution of the eigenvalues of the 
reduced density matrix. That is, the more uniformly dis- 
tributed the eigenvalues, the higher the entropy. In the 
PS phase, elements in the reduced density matrix related 
to the basis 1(3/3), which denotes two (3 atoms congregate 
together, are finite, while in the DW region, they are al- 
most zero. This fact leads to a larger entanglement with 
a block size larger than 2 in the PS phase, but otherwise 
in the DW phase. So the transition introduces a signifi- 
cant change into the value of the entanglement, and vice 
versa. From Fig. [T](c), we can see that E v {3) shares sim- 
ilar properties with E v (2). On the other hand, since the 
ground state is translational invariant, the entanglement 
satisfies the equation E v (l) = E V (L — I), Therefore, we 
have the same figures of E v (2) and E V (A) (Fig. [l](d)) for 
the 6-site system. 

Moreover, in the region I 6 [0, L/2] the entanglement 
is a non-decreasing function of I, as is shown in Fig. [5] 
for a 10-site system with N a = Np — 4 and U = 200. 
Therefore its first derivative develops a minimum at the 
critical point, as we can see from Fig. [21 Moreover, as the 
block size increases, the minimum point becomes sharper 
and sharper, exhibits a scaling behavior as shown in the 
inset of Fig. ® i.e., dE v (l) / dtp oc —I around the critical 
tp. 



IV. CHARGE ORDER PARAMETER AND 
PHASE DIAGRAM 

Though the entanglement can give us useful informa- 
tion about the phase diagram, the dominating configura- 
tions in different phases remains unknown. So it is impor- 
tant to study the structure factor in competing phases. 
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FIG. 3: (color online) The structure factor of DW as a func- 
tion of tp and various modes, i.e. quantized momentum. Here 
L = 10, N a = Np = 4, U = 200. 
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FIG. 5: (color online) The ground phase diagram in the U — tp 
plane in large U case. The dominating configuration in the 
up-region is DW and below the boundary line it is PS. Here 
L — 10, N a = Np. 
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FIG. 4: (color online) The structure factor of DW as a func- 
tion of tp and various modes, i.e. quantized momentum. Here 
L = 10, N a = N = 4, U = 50. 



Taking into account that the dominating configuration 
of (3 atoms is quite different in two phases, we introduce 
the following structure factor of DW of (3 atoms 



ScDw(?) = ^[ e " (H) (W 



{np) 



(15) 



where q = 2nir/L, n = 0, 1, • • • , L. In Fig. [3l we show 
that the structure factor as a function of tp for different 
modes for a system with L — 10 and a relatively large 
U = 200. The figure shows an obvious competition be- 
tween the two modes. In the small tp limit, i.e., when (3 
atom has very heavy mass, Scdw(<7 = 2tt/L) dominates, 
which indicates phase separation in this region [2(j. A 
careful scrutiny of the ground-state wavefunction finds 
that the configuration 

1/3,^,^,^,0,0,0,0,0,0) 



of P atoms is dominant. It is not difficult to interpret this 
result. In the small tp limit, the major contribution to 
the ground-state energy comes from the a atoms. So in 
order to have a lower energy, they need more free space 
because the energy of particles inside the Fermi surface 
is — 2cos(/c), k oc 1/L. Then f3 atoms will be pushed by a 
atoms to form clusters and phase separation occurs. For 
a finite size system, the translational symmetry is still 
preserved, while in the thermodynamic limit, especially 
in high dimensions, the symmetry in the ground state 
might be broken due to the very high potential energy 
between different configurations of (3 atoms. Then the 
system will be separated into two distinct regions macro- 
scopically. One is the solid-like region of (3 atoms, while 
the other is the liquid-like region of a atoms. The latter 
can be described by a model of iV Q atoms trapped in an 
infinite potential well with length L — Np + 1. The ground 
state is insulating and the energy is simply 



E 



-2 cos 



Jit 



L — Nr 



1 



(16) 



While if tp — > 1, Scdw(5 = Nir/L) exceeds <Scdw(<Z = 
2tt/L), which implies that (3 atoms distribute uniformly 
on the optical lattice. Then together with a atoms, the 
ground state becomes the so called DW state, which can 
be regarded as a solution of a and (3 atoms, as shown by 
the configuration 

\{3,a,o,f3,a,f3,o,a,f3, a). 

of a and (3 atoms. In its limiting case tp = 1, the model 
goes back to the traditional Hubbard model whose exci- 
tation spectrum is gapless, so the system is a conductor 
away from half- filling [13]. Therefore, different configura- 
tions dominate in different regions and the competition 



between them leads to a critical phenomenon. Accord- 
ing to this criterion, we can use the intersection of the 
structure factor of two modes to determine the transition 
point on the U — tp plane for a finite system. We plot 
the phase diagram on the U — tp plane in Fig. [5] for a 
10-site system with different filling N a = Np = 4, 6, 8. 

However, the results for a finite system is rather quali- 
tative. In order to have quantitative results for a real sys- 
tem, scaling analysis is crucial. For this purpose, we first 
estimate the scaling behavior of the ground-state energy 
in the critical region by the ED and DMRG method. [301] 
In Fig. [5J we show the scaling behavior of the ground- 
state energy at a given density n = 2/3. Results are 
obtained for systems with open boundary conditions via 
the DMRG method in which up to 150 states are kept in 
the finite algorithm. It is evident that the limiting energy 
is approached linearly with 1/L. A relation of the form 
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FIG. 6: The scaling behavior of the ground- state energy for 
two points near to the both sides of the critical point. Here 
n = 2/3. 



E Q (N) = E(po) +a/L, (17) 

where a is a constant, holds quite accurately in a large 
U region. Since the quantum critical phenomena is re- 
lated to the singularity in ground-state energy, the 1/L 
correlation in Eq. (fTTJ) actually implies that the phase 
boundary bears a similar scaling behavior. Based on this 
consideration, we take n = 2/3 as an example to show 
the scaling behavior of the phase boundary for both open 
and anti-periodic boundary conditions in Fig. \7\ 

It has been shown that, based on the variation princi- 
ples, one can obtain the lower and upper bounds of the 
phase boundary with different boundary conditions, such 
as periodic, anti-periodic, and open boundaries. From 
Fig. \7\ it is clearly shown that data with APBC give a 
lower bound while data with OBC give an upper bound 
on the transition point. Moreover, the extrapolated data 
based on the 1/L scaling of the two approaches for an 
infinite system agree with each other. This phenomenon 
is consistent with the fact that the physics in a real sys- 
tem should be independent of the boundary conditions. 
Moreover, we can also estimate errors in our extrapola- 
tion. We presented a final phase diagram with error bars 
smaller than the size of the symbols in Fig. [H 

In the small U region, we can see that there is a critical 
U on the [/-axis for the density n = 2/3. However, if n is 
reduced, the critical U tends to zero. In the low density 
limit, the phase boundary scales like tp oc U 2 in the small 
U region, which agrees with the results obtained by the 
Bosonization method [HJ excellently. 

In the large U region, the critical tp increases as U 
increases. We take a system of L — 12, N a = Np = 4 
as an example, and show the 1/U behavior of the phase 
boundary in Fig. [9] From the figure, we can see that 
the critical tp is proportional to 1/U in the large U limit. 
Moreover, Fig. [5] manifests that U will be saturated in 
the infinite U limit. That is, for a given concentration, 
there exists a saturation tp above which the phase sepa- 
ration will never happen regardless how large the on-site 
U is. Based on these physical intuition, the boundary 




FIG. 7: The scaling behavior of the phase boundary for both 
open boundary conditions (triangular lines) and anti-periodic 
boundary conditions (square lines). Here n = 2/3. 



line satisfies the relation 

tp = t s p + C/U, (18) 

where C is a constant, and both t% and C depend on the 
filling conditions. 

V. SINGLE-HOLE PROBLEM 

In this section, we give a rigorous proof that even for 
the case of one hole doping away from half-filling, the DW 
state is unstable to the PS state in the infinite U limit. 
If U is very large, the critical point then is approached 
linearly with 1/U. 

We first consider an odd-site sample with L — 2N a + 
1, N a — Np and infinite U. The space of DW is spanned 
by 2L basis: 

|ej) = |ai,/? 2 , ■•• °i ■■ .,a L -i,/3 L ),i 6 [1,L] 
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FIG. 8: The ground-state phase diagram of the AHM for two 
filling condtions: n = 2/3,2/5. 
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FIG. 9: 1/U behavior of the phase boundary in the large U 
region, as exemplified by a system with L = 12, N a = Np = 4. 



and 



|e t > = \fix, a 2) ■ ■ ■ o t _ L . . .,0 L -x,a L ),i G [L + 1,2L], 



where o, denotes a hole at site i. Then the Hamiltonian 
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FIG. 10: The periodical square potential wells used to de- 
scribe the dynamics of a hole in demixed phase. 
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(19) 



It can be solved exactly. The whole energy spectra of the 
system are given by 



E+=±Jl 



tp + 2tp cos(fcj), 



(20) 



k j = 3*/ L , J =0,1,. 



1 



The ground-state wavefunction for arbitrary odd L is 

1 



*n = 



2L 



2L f-f 

7 = 1 



(21) 



with the eigenenergy £dw = — 1 — tp. 

However, in the case of completely demixed phase, the 
whole space is spanned by L(L — 1) basis. The Hamilto- 
nian then describes the problem of single particle motion 
in periodic square potential wells (see Fig. fTU)) with dif- 
ferent hoping integrals in different region. Precisely, the 
Hamiltonian becomes 



H 



■^2d]d j+s , 

3,6 



(22) 



where dj and dj are hole creation and annihilation oper- 
ators in each potential well, and 



H = -tp d\dj + s 



(23) 



elsewhere. The ground state of this Hamiltonian in the 
thermodynamic limit is identical to the ground state of 



H 



-2+p 2 , 



< x < N a ; 
N a < x < Na, 



-It a + tap 
with periodic boundary conditions 

= *(.t + N a + N a ) 



(24) 



(25) 



For the latter, a bound state always exists for arbitrary 
well depth 2 — 2tg. Then if N a = Ng — ► oo, the ground- 
state energy is simply -2, which is obviously smaller than 
J?dw = — 1 — tp. For a finite system except L = 3, it can 
also be shown that Ep$(L) < Ejyw(L). For example if 
L = 5, N a = Ng — 2, we have 



E PS (L = 5) 



2t% 



(26) 



Therefore, in the infinite U limit, the DW state is unsta- 
ble to the PS state. Such a rigorous result is also valid 
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for a system of two species of hard-core bosonic atoms 
and Bose-Fermi mixtures with different hoping integrals 
in optical lattices. 

When the on-site U is very large but not infinite, the 
Hamiltonian |T]) can be approximated by [l3[ 



j=l <5=±1 a 
L 

+j£[S r S i+1 + *S«s£ +x 



/iTij-nj+i] ,(27) 



in which 

J = 



17 : 



<5 = 



2t a tp 



V = 



L a T 
8t a tp 



(28) 



Clearly if = tp, the above Hamiltonian becomes the 
< — J model. The ground state of the i — J model with a 
single-hole doping becomes the Nagaoka ferromagnetism 
3l| if the contribution from the kinetic energy in Eq. 
(|2T|) exceeds that from the spin-spin antiferromagnetic 
interaction. In order to study the condition of PS, we 
first suppose that the ground state of the system is phase 
separated. Then the ground-state energy can be approx- 
imated by 



Eps — —2 



(L-6)(1 + <5)J 



(29) 



On the other hand, if the ground state is in the DW state, 
the ground-state energy can be approximated by that of 
the XXZ chain. For the latter, the ground-state energy 
per bond has the form [33| 



VI. DISCUSSIONS 

Obviously, unlike the PS in the t — J model [15] and 
the extended Hubbard model [2g|, which is the conse- 
quence of attractive interaction between particles, the 
PS in the ground state of the AHM is driven by kinetic 
energy. Therefore, though our results are based on a one- 
dimensional model, the underlying physics is quite gen- 
eral for systems in any dimension. That is, in the large U 
limit, the dynamics of a system of two species of atoms 
at zero temperature is dominated by the light atoms. In 
order to have a lower energy, they need more free space 
to move. This mechanism forces heavy atoms to congre- 
gate together, so the latter becomes a solid-like object. 
In experiment, two separated regions are expected to be 
witnessed macroscopically. However, when tp — > t a , the 
dynamics of heavy atoms is comparable to that of light 
atoms, the exchange interaction drives the system into 
a DW state. Therefore, if we consider the PS state as 
a classical phase containing solid-like order and the DW 
state as a quantum region with liquid properties in the 
whole system, the transition reported in our work is just 
an example of a crossover from the classical region to the 
quantum region. 

Such an interesting transition is expected to be ob- 
served in the on-going experiments on optical lattices. 
We take a system consisting of two species of atoms (such 
as 6 Li (a) and 40 K (j3) with mp/m a ~ 20/3 ) as an exam- 
ple. Since the typical scattering length for alkaline atoms 
ranges from 40 to 100 aeohr [34|,|35j], and laser wavelength 
A = 852 nm ||. Then from Eq. {5]) we roughly estimate 
that PS phase can be observed when v > 0.4 according 
the phase diagram in Fig. [5J 



1 



exxz — 



— sinh i 



1 + 2 V 

2 ^—^ g2n0 

~ n=l 



1 



1 



(30) 



where cosh0 =1 + 5. Then the ground-stat energy of 
the DW phase becomes 



-1 - tp + (L — 2) Jexxz- 



(31) 



Here, in both £bw and -Bps, the finite-size correction to 
the ground-state energy is not taken into account, so the 
critical value is estimated approximately. Despite of this, 
the qualitative behavior of critical point is clear, i.e. 



tp ~ 1 — JL sinhi 



oo 



2 Z^ e 2n- 

n=l 



1 



1 



0(l/C/),(32) 



which means only if U 3> L, the DW state is unstable to 
PS state. For a given L, the phase boundary scales like 



tp = l- C/U, 

in the large U limit. Clearly, Eq. 
our previous result Eq. (|18p . 



(33) 

is consistent with 



VII. SUMMARY 

In summary, we have investigated the ground-state 
phase diagram of two species of fermionic atoms trapped 
in one-dimensional optical lattice. By using the ED 
method, we computed the block-block entanglement be- 
tween a local block and rest part for a small system. 
We obtained an intuitive picture of phase diagram of the 
ground state and found that the entanglement in the PS 
region is in general larger than that in the DW region 
for a finite system. Its first derivative develops a sharp 
downward peak and shows scaling behavior at the critical 
point. We also analyzed the structure factor of the DW of 
(3 atoms by the ED and DMRG method, and found that 
the competition between two different configurations in 
the ground-state wavefunction leads to a phase transition 
at the critical point. The global phase diagram was ob- 
tained from the careful scaling analysis for various-size 
systems and different boundary conditions. Therefore, 
we results firstly gave a quantitative description of the 
ground-state phase transition of the AHM away from the 
half-filling. Furthermore, we gave a rigorous proof that 
even for the case of a single hole doping, the DW state is 
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unstable to the PS in the infinite U limit. Such a rigor- 
ous conclusion clarifies the physical picture of the phase 
separation. 
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